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t^. ■ 1. Introduction. There are several equivalent formulations of the problem of quantizing an 

interacting field theory. The list includes canonical quantization, path integral (or functional) 
I techniques, stochastic quantization, "unified" methods such as the Batalin-Vilkovisky formalism, 

and techniques based on the realizations of field theories as low energy limits of string theory. 
The problem of obtaining an exact nonperturbative description of a given quantum field theory 
is most often a very difficult one. Perturbative techniques, on the other hand, are abundant, 
and common to all of the quantization methods mentioned above is that they admit particle 
ly-j ■ interpretations in this formalism. 

o : 

The basic physical quantities that one wishes to calculate in a relativistic (d+ l)-dimensional 
» , quantum field theory are the S-matrix elements 

Oh' 

s ba = ont (Mt)\Mt)) in (i) 

J> ! between in and out states at large positive time t. The scattering operator S is then defined by 

writing Q in terms of initial free particle (descriptor) states as 

. s ha =:(Mo)\s\Mo)) ■ (2) 

Suppose that the hamiltonian of the given field theory can be written as H = Hq + H', where 
Hq is the free part and H' the interaction hamiltonian. The time evolution of the in and out 
states are governed by the total hamiltonian H. They can be expressed in terms of descriptor 
states which evolve in time with Hq in the interaction picture and correspond to free particle 
states. This leads to the Dyson formula 



S = Texp^-i J tiAHj(t)j , (3) 

where T denotes time-ordering and Hj(t) = f d d x rii n t(x,t) is the interaction hamiltonian in 
the interaction picture, with fii n t(x,t) the interaction hamiltonian density, which deals with 
essentially free fields. This formula expresses S in terms of interaction picture operators acting 
on free particle states in (J2J) and is the first step towards Feynman perturbation theory. 

1 To appear in the Quantum Field Theory section of the Encyclopedia of Mathematical Physics, Elsevier, 2006. 
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For many analytic investigations, such as those which arise in renormalization theory, one 
is interested instead in the Green's functions of the quantum field theory, which measure the 
response of the system to an external perturbation. For definiteness, let us consider a free real 
scalar field theory in {d + l)-dimensions with lagrangian density 

C = \ <V> cHV - \ m 2 4? + £ int (4) 

where £; n t is the interaction lagrangian density which we assume has no derivative terms. The 
interaction hamiltonian density is then given by Ti m t = — Ant- Introducing a real scalar source 
J(x), we define the normalized "partition function" through the vacuum expectation values 



Z[J] = V, ' J (5) 





S[J] 


|o> 


<o 


|S[0]| 





where |0) is the normalized perturbative vacuum state of the quantum field theory given by 
(defined to be destroyed by all field annihilation operators), and 



S[J] = T exp J d d+1 x (£ int + J(x) cj){x)) J 



(6) 



from the Dyson formula. This partition function is the generating functional for all Green's 
functions of the quantum field theory, which are obtained from (jSJ) by taking functional deriva- 
tives with respect to the source and then setting J(x) = 0. Explicitly, in a formal Taylor series 
expansion in J one has 

oo . n n ,. 

z[j] = E M / dd+lxi J ^ ° {n) ^ ■•• >*«) . ( 7 ) 

n=0 U - i=l J 

whose coefficients are the Green's functions 

(n) ._ (0|T[exp(i/d^x £ int ) 0(x 1 )--^(x n )]|Q) 

( l '"" nj -" <0|T exp(i/d^xAnt)|0> • ^ 



It is customary to work in momentum space by introducing the Fourier transforms 
J(k)= I d d+1 x e ik ' x J(x) , & n \k 1 ,...,k n )=f[ f d d+1 x, e'^GW^,..,^) , (9) 

i=l 

in terms of which the expansion (J7J) reads 

oo . n n „ jd+1 1, 

z [ j ]=E^n T^j(-h)G^( kl ,...,k n ). (io) 

n=0 i=l J V j 

The generating functional Q1U|) can be written as a sum of Feynman diagrams with source 
insertions. Diagrammatically, the Green's function is an infinite series of graphs which can be 
represented symbolically as 
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where the n external lines denote the source insertions of momenta fcj and the bubble denotes 
the sum over all Feynman diagrams constructed from the interaction vertices of £i n t- 

The Green's functions can also be used to describe scattering amplitudes, but there are two 
important differences between the graphs and those which appear in scattering theory. In 
the present case, external lines carry propagators, i.e. the free field Green's functions 

A(x-y) = <0|T [«/>(*) ^)] \0) = (x\(n +m y 1 \y) = /(j^ggr p2-^2 +ie e ~ iHX ~ y) 

(12) 

where e — ► + regulates the mass shell contributions, and their momenta ki are off-shell in 
general (kf 7^ to 2 ). By the LSZ theorem, the S-matrix element is then given by the multiple 
on-shell residue of the Green's function in momentum space as 



n I 




x G( n + m )(-k' 1 ,...,-k , n ,k 1 ,...,k l ) , (13) 

where icj, ic^ are the residues of the corresponding particle poles in the exact two-point Green's 
function. 

This article deals with the formal development and computation of perturbative scattering 
amplitudes in relativistic quantum field theory, along the lines outlined above. Initially we deal 
only with real scalar field theories of the sort © in order to illustrate the concepts and technical 
tools in as simple and concise a fashion as possible. These techniques are common to most 
quantum field theories. Fermions and gauge theories are then separately treated afterwards, 
focusing on the methods which are particular to them. 



2. Diagrammatics. The pinnacle of perturbation theory is the technique of Feynman 
diagrams. Here we develop the basic machinery in a quite general setting and use it to analyse 
some generic features of the terms comprising the perturbation series. 

Wick's theorem. The Green's functions (JSJ) are defined in terms of vacuum expectation values 
of time-ordered products of the scalar field 4>(x) at different spacetime points. Wick's theorem 
expresses such products in terms of normal-ordered products, defined by placing each field 
creation operator to the right of each field annihilation operator, and in terms of two-point 
Green's functions ((12)) of the free field theory (propagators). The consequence of this theorem 
is the haffnian formula 



<O|T[0(a;i)-..^(x n )] |o) 







e n {^w^-i))^^))} |o) 



n 



2k - 1 , 



n 



2k 



(14) 



The formal Taylor series expansion of the scattering operator S may now be succinctly sum- 
marized into a diagrammatic notation by using Wick's theorem. For each spacetime integration 
J d d+l Xi we introduce a vertex with label i, and from each vertex there emanates some lines 
corresponding to field insertions at the point Xj. If the operators represented by two lines ap- 
pear in a two-point function according to (|14|). i.e. they are contracted, then these two lines are 
connected together. The S operator is then represented as a sum over all such Wick diagrams, 
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bearing in mind that topologically equivalent diagrams correspond to the same term in S. Two 
diagrams are said to have the same pattern if they differ only by a permutation of their vertices. 
For any diagram D with n(D) vertices, the number of ways of interchanging vertices in n(D)\. 
The number of diagrams per pattern is always less than this number. The symmetry number 
S(D) of D is the number of permutations of vertices that give the same diagram. The number 
of diagrams with the pattern of D is then n(D)\/S(D). 

In a given pattern, we write the contribution to S of a single diagram D as n /^j ° #(2)) ° , 
where the combinatorial factor comes from the Taylor expansion of S, the colons denote normal 
ordering of quantum operators, and ° 0(D) ° contains spacetime integrals over normal-ordered 
products of the fields. Then all diagrams with the pattern of D contribute ° 0(D) ° /S(D) to S. 
Only the connected diagrams D r , r S N (those in which every vertex is connected to every other 
vertex) contribute and we can write the scattering operator in a simple form which eliminates 
contributions from all disconnected diagrams as 

Feynman rules. Feynman diagrams in momentum space are defined from the Wick diagrams 
above by dropping the labels on vertices (and also the symmetry factors S(D) _1 ) : and by la- 
belling the external lines by the momenta of the initial and final particles that the corresponding 
field operators annihilate. In a spacetime interpretation, external lines represent on-shell phys- 
ical particles while internal lines of the graph represent off-shell virtual particles (k 2 ^ m 2 ). 
Physical particles interact via the exchange of virtual particles. An arbitrary diagram is then 
calculated via the Feynman rules 

f d d+1 p i 
J (2n) d + 1 p 2 -m 2 +ie ' 

ig (2vr) d+1 6( d+1 \ Pl + • • • + p n ) 

(16) 



Irreducible Green's functions. A one-particle irreducible (1PI) or proper Green's function is 
given by a sum of diagrams in which each diagram cannot be separated by cutting one internal 
line. In momentum space, it is defined without the overall momentum conservation delta- 
function factors and without propagators on external lines. For example, the two particle 1PI 
Green's function 

L =: E(Jfe) 

(17) 

is called the self-energy. If G(k) is the complete two-point function in momentum space, then 




for a monomial interaction C\ — — ^ n 



'int — „| 
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one has 



G(k) := ± 



k 



+ 




+ 



i 

k 2 - m 2 - ( 18 ) 
and thus it suffices to calculate only 1PI diagrams. 

The 1PI effective action, defined by the Legendre transformation r[^>] := — ilnZ[J] — 
/ d d+1 x J(x)<j>(x) of ©, is the generating functional for proper vertex functions and it can 
be represented as a functional of only the vacuum expectation value of the field <fi, i.e. its clas- 
sical value. In the semi-classical (WKB) approximation, the one-loop effective action is given 
by 

T[<j>] = S[<f>] + f Tr ln(l + A V" [<(>]) + O {h 2 ) 

00 (_-\\n n r 

= s[cp\ { -^- n / dd+i ^ a ^ - ^+1) ^'[^(^+1)] + ° ( h2 ) ( i9 ) 

n=l n i=l J 

where we have denoted S[<p] = J d d+l x C and V[<p] = —C mt , and for each term in the infinite 
series we define x n +i := x\. The first term in (|19|) is the classical contribution and it can be 
represented in terms of connected tree diagrams. The second term is the sum of contributions 
of one-loop diagrams constructed from n propagators — iA(x — y) and n vertices —iV" [<f>]. 
The expansion may be carried out to all orders in terms of connected Feynman diagrams, and 
the result of the above Legendre transformation is to select only the one-particle irreducible 
diagrams and to replace the classical value of <f> by an arbitrary argument. All information 
about the quantum field theory is encoded in this effective action. 

Parametric representation. Consider an arbitrary proper Feynman diagram D with n internal 
lines and v vertices. The number I of independent loops in the diagram is the number of 
independent internal momenta in D when conservation laws at each vertex have been taken 
into account, and it is given by £ = n + 1 — v . There is an independent momentum integration 
variable k{ for each loop, and a propagator for each internal line as in The contribution 

of D to a proper Green's function with r incoming external momenta pi, with X^i=iK = 0) i s 
given by 

'»w - H n / ^ ^tt- n s (d+l} m - <»> 

i=i 1 j=i 

where V(D) contains all contributions from the interaction vertices of £- m t, and Pj (resp. Kj) is 
the sum of incoming external momenta p\. (resp. internal momenta ki ) at vertex j with respect 
to a fixed chosen orientation of the lines of the graph. After resolving the delta- functions in 
terms of independent internal loop momenta k\ , . . . , kg and dropping the overall momentum 
conservation delta- function along with the symmetry and vertex factors in (|2U|). one is left with 
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a set of momentum space integrals 

where a,j(k,p) are functions of both the internal and external momenta. 

It is convenient to exponentiate propagators using the Schwinger parametrization 

dctj e iQ ^ +ie ) , (22) 



a,j + i e 



o 



and after some straightforward manipulations one may write the Feynman parametric formula 



n ^Tv t = <» - D' n J (^) 

7 — 1 .7 — 1 n 



J =1 . 

where D^,(k;a,p) := ^ • Oj[ay(fc,p) + ie] is generically a quadratic form 



D v (k;a,p) = - £ A* • QiMh + ^ L i(p) ' ki + X (P 2 ) ■ ( 24 ) 



The positive symmetric matrix is independent of the external momenta pi, invertible, and has 
non-zero eigenvalues Qi, ■ ■ ■ ,Qe- The vectors are linear combinations of the pj, while \(p 2 ) 
is a function of only the Lorentz invariants pf. After some further elementary manipulations, 
the loop diagram contribution (|21|) may be written as 



n f 1 1 f d d+1 k- 

5 Ei ^ 2 + a (p 2 ) - 1 Eij ^(p) • Q-HukHp))^ ■ (25) 



Finally, the integrals over the loop momenta ki may be performed by Wick rotating them to 
euclidean space and using the fact that the combination of £ integrations in M. d+1 has 0((d+ 1)£) 
rotational invariance. The contribution from the entire Feynman diagram T> thereby reduces to 
the calculation of the parametric integrals 



r(n-^) » } i 1 H!-E,«ii 

IM = (2n)^^ n J ^ II W ; —SIM 

(2vr) 2 1 ^ =1 o (a^-IE^^W-Q-H^Cp)) 2 

(26) 

where T(s) is the Euler gamma-function. 



3. Regularization. The parametric representation 1)26(1 is generically convergent when 
2n — (d + 1)£ > 0. When divergent, the infinities arise from the lower limits of integration 
a.j — > 0. This is just the parametric representation of the large k divergence of the original 
Feynman amplitude (|2U() . Such ultraviolet divergences plague the very meaning of a quantum 
field theory and must be dealt with in some way. We will now quickly tour the standard methods 
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of ultraviolet regularization for such loop integrals, which is prelude to the renormalization 
program that removes the divergences (in a renormalizable field theory). Here we consider 
regularization simply as a means of justification for the various formal manipulations that are 
used in arriving at expressions such as (|26[). 

Momentum cutoff. Cutoff regularization introduces a mass scale A into the quantum field 
theory and throws away the Fourier modes of the fields for spatial momenta k with \k\ > A. 
This regularization spoils Lorentz invariance. It is also non-local. For example, if we restrict to 
a hypercube in momentum space, so that \ki\ < A for i = 1, . . . , d, then Jj fc | <A &J?d e lk x = 

rii=i sin(Ax*)/7r x l which is a delta-function in the limit A — > oo but is non-local for A < oo. 
The regularized field theory is finite order by order in perturbation theory and depends on the 
cutoff A. 

Lattice regularization. We can replace the spatial continuum by a lattice £ of rank d and define 
a lagrangian on £ by 

Ac = ^E^ 2 + J E (27) 

ies(£) (i,j)eL(£) teS(£) 

where S(£) is the set of sites i of the lattice on each of which is situated a time-dependent 
function fa, and L(£) is the collection of links connecting pairs of nearest-neighbour sites 
i,j on £. The regularized field theory is now local, but still has broken Lorentz invariance. In 
particular, it suffers from broken rotational symmetry. If £ is hypercubic with lattice spacing 
a, i.e. £ = (Za) rf , then the momentum cutoff is at A = a" 1 . 

Pauli-Villars regularization. We can replace the propagator i(k 2 — m 2 + ie)" 1 by i (k 2 — 



rn 



+ ie) _1 + i J2j=i c j(k 2 ~ + ie) -1 , where the masses Mj 3> m are identified with the 
momentum cutoff as min{Mj} = A — > oo. The mass-dependent coefficients Cj are chosen to 
make the modified propagator decay rapidly as (/c 2 ) _Ar_1 at k — > oo, which gives the N equations 
(to 2 )* + Yjj c j( M jT = 0, i = 0,1,... ,N - 1. This regularization preserves Lorentz invariance 
(and other symmetries that the field theory may possess) and is local in the following sense. 
The modified propagator can be thought of as arising through the alteration of the lagrangian 
density (@J) by N additional scalar fields (fj of masses Mj with 

N 

£pv = \ 8^ - \ m 2 cf> 2 + ( \ &m d^Vj ~ \ M- + Ant m (28) 

where $ := 4>+Ylj \f^i fj- The contraction of the $ field thus produces the required propagator. 
However, the Cj's as computed above are generically negative numbers and so the lagrangian 
density (|28j) is not hermitian (as $ It is possible to make (|28[) formally hermitian by 

redefining the inner product on the Hilbert space of physical states, but this produces negative 
norm states. This is no problem at energy scales E <C Mj on which the extra particles decouple 
and the negative probability states are invisible. 

Dimensional regularization. Consider a euclidean space integral f d 4 k (k 2 + a 2 )~ r arising after 
Wick rotation from some loop diagram in (3 + l)-dimensional scalar field theory. We replace 
this integral by its D-dimensional version 



d ° k _ ^(- 2 )° 2 ' r(r m (99) 
{k 2 +a 2 ) r (r-l)\ { 2) ■ [ > 
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This integral is absolutely convergent for D < 2r. We can analytically continue the result of 
this integration to the complex plane D € C. As an analytic function, the only singularities of 
the Euler function T(z) are poles at z = 0, —1, —2, .... In particular, T(z) has a simple pole at 
z = of residue 1. If we write D = 4 + e with |e| — > 0, then the integral (|29j) is proportional to 
r(r — 2 — |) and e plays the role of the regulator here. This regularization is Lorentz invariant 
(in D dimensions) and is distinguished as having a dimensionless regularization parameter e. 
This parameter is related to the momentum cutoff A by e _1 = ln(A/m), so that the limit e — ► 
corresponds to A — > oo. 

Infrared divergences. Thus far we have only considered the ultraviolet behaviour of loop 
amplitudes in quantum field theory. When dealing with massless particles (m = in @) one 
has to further worry about divergences arising from the k — > regions of Feynman integrals. 
After Wick rotation to euclidean momenta, one can show that no singularities arise in a given 
Feynman diagram as some of its internal masses vanish provided that all vertices have superficial 
degree of divergence d + 1, the external momenta are not exceptional (i.e. no partial sum of 
the incoming momenta pi vanishes), and there is at most one soft external momentum. This 
result assumes that renormalization has been carried out at some fixed euclidean point. The 
extension of this property when the external momenta are continued to physical on-shell values 
is difficult. The Kinoshita-Lee-Nauenberg theorem states that, as a consequence of unitarity, 
transition probabilities in a theory involving massless particles are finite when the sum over all 
degenerate states (initial and final) is taken. This is true order by order in perturbation theory 
in bare quantities or if minimal subtraction renormalization is used (to avoid infrared or mass 
singularities in the renormalization constants). 



4. Fermion Fields. We will now leave the generalities of our pure scalar field theory and start 
considering the extensions of our previous considerations to other types of particles. Henceforth 
we will primarily deal with the case of (3 + l)-dimensional spacetime. We begin by indicating 
how the rudiments of perturbation theory above apply to the case of Dirac fermion fields. The 
lagrangian density is 

C F = ^(i$ -m)i) + £ (30) 

where tp are four-component Dirac fermion fields in (3 + l)-dimensions, ip := ^7° and ^ = 7^ 
with 7^ the generators of the Clifford algebra {7^, 7^} = 2n^ v . The lagrangian density C 
contains couplings of the Dirac fields to other field theories, such as the scalar field theories 
considered previously. 

Wick's theorem for anticommuting Fermi fields leads to the pfaffian formula 

, n = 2k - 1 , 



(0|T[V(1)---V(n)] |0) 



hi E sgn(vr) ]1 <0|T[V(7r(2i-l))^(vr(2i))] |0> , n = 2k 



2 k k\ 

TT€S 2 k « = 1 



(31) 

where for compactness we have written in the argument of ip(i) the spacetime coordinate, the 
Dirac index, and a discrete index which distinguishes ip from %p. The non- vanishing contractions 
in (|31|) are determined by the free fermion propagator 

A F (x-y) = <0|T[V»(x)?(y)] |0> =(x\(i$ -m)-'\y) = i j ^ ^^2% ■ 

(32) 

Perturbation theory now proceeds exactly as before. Suppose that the coupling lagrangian 
density in (|3U|) is of the form C! = ip(x) V(x) ijj(x). Both the Dyson formula © and the 
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diagrammatic formula (|15|) are formally the same in this instance. For example, in the for- 
mal expansion in powers of J d 4 x £ , the vacuum-to- vacuum amplitude (the denominator in 
((5J)) will contain field products of the form YYi=i I (0| r ^['4 , ( x i) V( x %) VK 3 ^)] |0) which cor- 
respond to fermion loops. Before applying Wick's theorem the fields must be rearranged as 
tr niLi V( x i) VK 3 ^) i)(. x i+i) (with x n+ \ := xi), where tr is the 4x4 trace over spinor indices. 
This reordering introduces the familiar minus sign for a closed fermion loop, and one has 




Feynman rules are now described as follows. Fermion lines are oriented to distinguish a 
particle from its corresponding antiparticle, and carry both a four-momentum label p as well as 
a spin polarization index r = 1,2. Incoming fermions (resp. antifermions) are described by the 
wavefunctions u p r ^ (resp. Vp V ^), while outgoing fermions (resp. antifermions) are described by 
the wavefunctions Up^ (resp. i>p ). Here and v p r ^ are the classical spinors, i.e. the positive 
and negative energy solutions of the Dirac equation (j$ — m)u p r ^ = (p 1 + m)v p r ^ = 0. Matrices 
are multiplied along a Fermi line, with the head of the arrow on the left. Closed fermion loops 
produce an overall minus sign as in (|33[) . and the multiplication rule gives the trace of Dirac 
matrices along the lines of the loop. Unpolarized scattering amplitudes are summed over the 
spins of final particles and averaged over the spins of initial particles using the completeness 
relations for spinors 

^ff=^+ m , J2vPvM=t-m, (34) 
r=l,2 r=l,2 

leading to basis independent results. Polarized amplitudes are computed using the spinor bilin- 
ears Up ^ p u p — v p ^ p v p = 2p p o , Up u p = —Vp v p = 2m o , and Up v p = 0. 

When calculating fermion loop integrals using dimensional regularization, one utilizes the 
Dirac algebra in D dimensions 

D , 

{2-D)j> , 
4p • k + (D - 4) p 1 $ , 
-2^j*-(Z>-4)j*^ , 
, trW = 4rf , 

4 (rf rf - rf if + rf ia rf) . (35) 

Specific to D = 4 dimensions are the trace identities 

tr 7 5 = tr 7^ 7^ 7 5 = , tr j p Y Y 7 CT 7 5 = -4 i e pupa (36) 

where 7 5 := i7°7 1 7 2 7 3 . Finally, loop diagrams evaluated with the fermion propagator (|32|) 
require a generalization of the momentum space integral ()29j) given by 

r d D k i 

J Wf 5 {ki + 2k.p + a?+\e) r " (2vr) B (r - 1)! ( a 2 _ p 2 + ie y-# ' 



ini> 7 M = 

7^#7m = 
trl = 4 , tr 7 Ml •••7 M2fe - 1 = 



i(-.)fF(r-f) 
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From this formula we can extract expressions for more complicated Feynman integrals which are 
tensorial, i.e. which contain products of momentum components k^ in the numerators of their 
integrands, by differentiating (|37|) with respect to the external momentum p^. 



5. Gauge Fields. The issues we have dealt with thus far have interesting difficulties when 
dealing with gauge fields. We will now discuss some general aspects of the perturbation expansion 
of gauge theories using as prototypical examples quantum electrodynamics (QED) and quantum 
chromodynamics (QCD) in four spacetime dimensions. 

Quantum electrodynamics. Consider the QED lagrangian density 

Cq^n = -\F llv F^ + \ix 2 A ll A> l + ^{i$ -e4 -m)i/> (38) 

where is a U(l) gauge field in (3 + l)-dimensions and F^ u = d^A u — d u A^ is its field 
strength tensor. We have added a small mass term fj, — > for the gauge field which at the 
end of calculations should be taken to vanish in order to describe real photons (as opposed 
to the soft photons described by (|38|)1. This is done in order to cure the infrared divergences 
generated in scattering amplitudes due to the masslessness of the photon, i.e. the long range 
nature of the electromagnetic interaction. The Bloch-Nordsieck theorem in QED states that 
infrared divergences cancel for physical processes, i.e. for processes with an arbitrary number of 
undetectable soft photons. 

Perturbation theory proceeds in the usual way via the Dyson formula, Wick's theorem, and 
Feynman diagrams. The gauge field propagator is given by 

(OIT^A^)] |0> = (,|[^ C + M 2 ) -dM~ X \y) = ij ^ JZfJle 

(39) 

and is represented by a wavy line. The fermion-fermion-photon vertex is 

< = - ie7^ • 



(40) 



An incoming (resp. outgoing) soft photon of momentum k and polarization r is described by 
the wavefunction e^\k) (resp. e^\k)*), where the polarization vectors (k), r = 1,2,3 
solve the vector field wave equation (□ + [i 2 )A^ = d^A 11 = and obey the orthonormality and 
completeness conditions 

e« (k)* ■ e« (k) = —S rs , 4 ] ( k ) ^ (*0* = + ^T L , (41) 

r=l ^ 

along with k ■ e^'\k) = 0. All vector indices are contracted along the lines of the Feynman 
graph. All other Feynman rules are as previously. 

Quantum chromodynamics. Consider nonabelian gauge theory in (3 + l)-dimensions minimally 
coupled to a set of fermion fields ip , A = l,...,Nf each transforming in the fundamental 
representation of the gauge group G whose generators T a satisfy the commutation relations 
[T a ,T b ] = f abc T c . The lagrangian density is given by 

Nf 

£qcd = -\ F; v F a ^ + ± (d,A a n 2 + dtfLPri + J2^ A (ip- m A )^ A , (42) 

A=l 
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where F^ v 



d v A% + f abc AlAl and D„ 



df, + ieR(T a )A^ with R the pertinent 
representation of G (R(T a )b c = f a bc for the adjoint representation and R(T a ) = T a for the 
fundamental representation). The first term is the Yang-Mills lagrangian density, the second 
term is the covariant gauge-fixing term, and the third term contains the Faddeev-Popov ghost 
fields rj which transform in the adjoint representation of the gauge group. 

Feynman rules are straightforward to write down and are given by 



b.u 



a,iM 



-iS' 



ab 



Vtiu + (a - 1) 
k 2 + ie 



c,A,fc 




e f abc [v\n(Pv ~ K) + V^ilx ~ Px) + rjv\(kp - q»)] 



- i e 2 [f eab / ecd (^ A r/ vp - ripp r} vX ) 

+f eac f ebd ^ ^ _ ^ ^ +f ead febc^ ^ _ ^ ^ 



i<5 afe 
fc 2 + ie 

(43) 

where wavy lines represent gluons and dashed lines represent ghosts. Feynman rules for the 
fermions are exactly as before, except that now the vertex (j40j) is multiplied by the colour 
matrix T a . All colour indices are contracted along the lines of the Feynman graph. Colour 
factors may be simplified by using the identities 

Tr R a R b = di^B C2 (jq 6 ab ^ R a R a = ^ R a R b R a = ^{R) ~ \ C 2 (G)) R b (44) 

where R a := R(T a ) and C2(R) is the quadratic Casimir invariant of the representation R 
(with value C2(G) in the adjoint representation). For G = SU(N), one has C2{G) = N and 
C2(N) = (N 2 — 1)/2N for the fundamental representation. 

The cancellation of infrared divergences in loop amplitudes of QCD is far more delicate than 
in QED, as there is no analog of the Bloch-Nordsieck theorem in this case. The Kinoshita-Lee- 
Nauenberg theorem guarantees that, at the end of any perturbative calculation, these divergences 
must cancel for any appropriately defined physical quantity. However, at a given order of pertur- 
bation theory, a physical quantity typically involves both virtual and real emission contributions 
that are separately infrared divergent. Already at two- loop level these divergences have a highly 
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intricate structure. Their precise form is specified by the Catani colour-space factorization for- 
mula which also provides an efficient way of organising amplitudes into divergent parts, which 
ultimately drop out of physical quantities, and finite contributions. 

The computation of multi-gluon amplitudes in nonabelian gauge theory is rather complicated 
when one uses polarization states of vector bosons. A much more efficient representation of 
amplitudes is provided by adopting a helicity (or circular polarization) basis for external gluons. 
In the spinor-helicity formalism, one expresses positive and negative helicity polarization vectors 
in terms of massless Weyl spinors \k^} := s(l± 75)1^ = b(1± 75)^ through 

where q is an arbitary null reference momentum which drops out of the final gauge invariant 
amplitudes. The spinor products are crossing symmetric, antisymmetric in their arguments, and 
satisfy the identities 

(h I k/) (kf I k i ) = 2ki ■ kj , 

(K I ( k T I k r) = (K I kt) (kj- I fc+) + (K I k+)(kj I k+) . (46) 

Any amplitude with massless external fermions and vector bosons can be expressed in terms 
of spinor products. Conversely, the spinor products offer the most compact representation of 
helicity amplitudes which can be related to more conventional amplitudes described in terms of 
Lorentz invariants. For loop amplitudes, one uses a dimensional regularization scheme in which 
all helicity states are kept four-dimensional and only internal loop momenta are continued to 
D = 4 + e dimensions. 



6. Computing Loop Integrals. At the very heart of perturbative quantum field theory is 
the problem of computing Feynman integrals for multi-loop scattering amplitudes. The inte- 
grations typically involve serious technical challenges and for the most part are intractable by 
straightforward analytical means. We will now survey some of the computational techniques that 
have been developed for calculating quantum loop amplitudes which arise in the field theories 
considered previously. 

Asymptotic expansion. In many physical instances one is interested in scattering amplitudes in 
certain kinematical limits. In this case one may perform an asymptotic expansion of multi-loop 
diagrams whose coefficients are typically non-analytic functions of the perturbative expansion 
parameter H. The main simplification which arises comes from the fact that the expansions 
are done before any momentum integrals are evaluated. In the limits of interest, Taylor series 
expansions in different selected regions of each loop momentum can be interpreted in terms of 
subgraphs and co-subgraphs of the original Feynman diagram. 

Consider a Feynman diagram 2) which depends on a collection {Q{\ of large momenta (or 
masses), and a collection {mi,qi} of small masses and momenta. The prescription for the large 
momentum asymptotic expansion of D may be summarized in the diagrammatic formula 

lim D(Q;m,q) = ^ {D/d)(m,q) -k (T{ 

mo, go }*>) (Q^tbQo) j (47) 

DCS 

where the sum runs through all subgraphs 5 of ID which contain all vertices where a large 
momentum enters or leaves the graph and is one-particle irreducible after identifying these 
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vertices. The operator Ttm^^x performs a Taylor series expansion before any integration is 
carried out, and the notation (2)/fl) * (T^ mo qj) yd) indicates that the subgraph d C D is replaced 
by its Taylor expansion in all masses and external momenta of D that do not belong to the set 
{Qi}. The external momenta of d which become loop momenta in D are also considered to 
be small. The loop integrations are then performed only after all these expansions have been 
carried out. The diagrams D/d are called co-subgraphs. 

The subgraphs become massless integrals in which the scales are set by the large momenta. 
For instance, in the simplest case of a single large momentum Q one is left with integrals 
over propagators. The co-subgraphs may contain small external momenta and masses, but the 
resulting integrals are typically much simpler than the original one. A similar formula is true 
for large mass expansions, with the vertex conditions on subdiagrams replace by propagator 
conditions. For example, consider the asymptotic expansion of the two-loop double bubble 
diagram 





O 




(48) 

in the region q 2 <C m 2 , where m is the mass of the inner loop. The subgraphs (to the right of 
the stars) are expanded in all external momenta including q and reinserted into the fat vertices 
of the co-subgraphs (to the left of the stars). Once such asymptotic expansions are carried out, 
one may attempt to reconstruct as much information as possible about the given scattering 
amplitude by using the method of Pade approximation which requires knowledge of only part of 
the expansion of the diagram. By construction, the Pade approximation has the same analytic 
properties as the exact amplitude. 

Brown- Feynman reduction. When considering loop diagrams which involve fermions or gauge 
bosons, one encounters tensorial Feynman integrals. When these involve more than three distinct 
denominator factors (propagators), they require more than two Feynman parameters for their 
evaluation and become increasingly complicated. The Brown-Feynman method simplifies such 
higher-rank integrals and effectively reduces them to scalar integrals which typically require 
fewer Feynman parameters for their evaluation. 

To illustrate the idea behind this method, consider the one-loop rank 3 tensor Feynman 
integral 

tuuX f d D k k»k u k x 

J (2ir) D k 2 {k 2 - fi 2 )(q-k) 2 ((k-q) 2 + fi 2 )(k 2 + 2k-p) V ; 

where p and q are external momenta with the mass-shell conditions p 2 = (p — q) 2 = m 2 . By 
Lorentz invariance, the general structure of the integral (|49f) will be of the form 

where a} w , b^ v are tensor- valued functions and &" a vector- valued function of p and q. The 
symmetric tensor is chosen to project out components of vectors transverse to both p and 
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q, i.e. p^ s^ v = q^ s^ v = 0, with the normalization stf = D — 2. Solving these constraints leads 
to the explicit form 



m 2 qK q v + q 2 p» p v - (p ■ q) (g^ p v + p» q v ) 
m A q* — [p ■ qy 



To determine the as yet unknown functions a^ u , b^ v and &" above, we first contract both 
sides of the decomposition (|5U|) with p^ and q^ to get 

2p x J^ x = 2m 2 a^ + 2(p-q)b^ , 2q x J^ A = 2(p • q) oT + 2g 2 \f v . (52) 

Inside the integrand of (|49|). we then use the trivial identities 

2k - p = (k 2 + 2k ■ p) - k 2 , 2q ■ k = k 2 + q 2 - (k - q) 2 (53) 

to write the left-hand sides of (|52j) as the sum of rank 2 Feynman integrals which, with the 
exception of the one multiplied by q 2 from (|53|) . have one less denominator factor. This formally 
determines the coefficients a^ v and b^ v in terms of a set of rank 2 integrations. The vector 
function is then found from the contraction 

J^ v = Vv a^ v + q v b» u + (D - 2) d 1 . (54) 

This contraction eliminates the k 2 denominator factor in the integrand of (|49|) and produces a 
vector-valued integral. Solving the system of algebraic equations (|52|) and ()54|) then formally 
determines the rank 3 Feynman integral (|49|) in terms of rank 1 and rank 2 Feynman integrals. 
The rank 2 Feynman integrals thus generated can then be evaluated in the same way by writing 
a decomposition for them analogous to 1)50(1 and solving for them in terms of vector- valued and 
scalar-valued Feynman integrals. Finally, the rank 1 integrations can be solved for in terms of a 
set of scalar-valued integrals, most of which have fewer denominator factors in their integrands. 

Generally, any one-loop amplitude can be reduced to a set of basic integrals by using the 
Passarino-Veltman reduction technique. For example, in supersymmetric amplitudes of gluons 
any tensor Feynman integral can be reduced to a set of scalar integrals, i.e. Feynman integrals 
in a scalar field theory with a massless particle circulating in the loop, with rational coefficients. 
In the case of J\f = 4 supersymmetric Yang-Mills theory, only scalar box integrals appear. 

Reduction to master integrals. While the Brown-Feynman and Passarino-Veltman reductions 
are well-suited for dealing with one-loop diagrams, they become rather cumbersome for higher 
loop computations. There are other more powerful methods for reducing general tensor integrals 
into a basis of known integrals called master integrals. Let us illustrate this technique on a scalar 
example. Any scalar massless two-loop Feynman integral can be brought into the form 

where Aj are massless scalar propagators depending on the loop momenta k, k' and the external 
momenta p±, . . . ,p n , and Sj are scalar products of a loop momentum with an external momentum 
or of the two loop momenta. The topology of the corresponding Feynman diagram is uniquely 
determined by specifying the set Ai, . . . , At of t distinct propagators in the graph, while the 
integral itself is specified by the powers lj > 1 of all propagators, by the selection Ei, . . . , S 9 of 
q scalar products and by their powers rij > 0. 
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The integrals in a class of diagrams of the same topology with the same denominator di- 
mension r = Ylj lj an d same total scalar product number s = Yli n i are related by various 
identities. One class follows from the fact that the integral over a total derivative with respect 
to any loop momentum vanishes in dimensional regularization as J jk^rb = 0' where J(k) 

is any tensorial combination of propagators, scalar products and loop momenta. The resulting 
relations are called integration by parts identities and for two-loop integrals can be cast into the 
form 

d D k f d D k> dfjkj^p) [ d D k [ d D k> dfjkj^p) 

~ v aUT. = = / j—r jr / t—^ — , (56) 



(2tt) d J (2n) D dk» J (2tt) d J (2n) D dk'» 

where f(k,k',p) is a scalar function containing propagators and scalar products, and v^ 1 is any 
internal or external momentum. For a graph with I loops and n independent external momenta, 
this results in a total of I (n + £) relations. 

In addition to these identities, one can also exploit the fact that all Feynman integrals 
1)55(1 are Lorentz scalars. Under an infinitesimal Lorentz transformation p M — > p^ + 5p^, with 
5p^ = p y SeJ*, SeJ 1 = —Se^, one has the invariance condition I(p + dp) = I(p), which leads to 
the linear homogeneous differential equations 



E 



*5b-rf5r)'M- a - (57> 



This equation can be contracted with all possible antisymmetric combinations of Pi^Pju to yield 
linearly independent Lorentz invariance identities for (|55|). 

Using these two sets of identities, one can either obtain a reduction of integrals of the type 
(|55[) to those corresponding to a small number of simpler diagrams of the same topology and 
diagrams of simpler topology (fewer denominator factors), or a complete reduction to diagrams 
with simpler topology. The remaining integrals of the topology under consideration are called 
irreducible master integrals. These momentum integrals cannot be further reduced and have to 
be computed by different techniques. For instance, one can apply a Mellin-Barnes transformation 
of all propagators given by 



1 oo 

1 1 f dz 



(k 2 + a) 1 (l-l) ] . J 2vri (k 2 ) l+z 



T(l + z)T(-z), (58) 



where the contour of integration is chosen to lie to the right of the poles of the Euler function 
r(/ + z) and to the left of the poles of T(—z) in the complex z-plane. Alternatively, one may 
apply the negative dimension method in which D is regarded as a negative integer in intermediate 
calculations and the problem of loop integration is replaced with that of handling infinite series. 
When combined with the above methods, it may be used to derive powerful recursion relations 
among scattering amplitudes. Both of these techniques rely on an explicit integration over the 
loop momenta of the graph, their differences occuring mainly in the representations used for the 
propagators. 

The procedure outlined above can also be used to reduce a tensor Feynman integral to scalar 
integrals, as in the Brown- Feynman and Passarino-Veltman reductions. The tensor integrals are 
expressed as linear combinations of scalar integrals of either higher dimension or with propaga- 
tors raised to higher powers. The projection onto a tensor basis takes the form (|55|) and can 
thus be reduced to master integrals. 
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7. String Theory Methods. The realizations of field theories as the low-energy limits of 
string theory provides a number of powerful tools for the calculation of multi-loop amplitudes. 
They may be used to provide sets of diagrammatic computational rules, and they also work well 
for calculations in quantum gravity. In this final part we shall briefly sketch the insights into 
perturbative quantum field theory that are provided by techniques borrowed from string theory. 

String theory representation. String theory provides an efficient compact representation of 
scattering amplitudes. At each loop order there is only a single closed string diagram, which 
includes within it all Feynman graphs along with the contributions of the infinite tower of massive 
string excitations. Schematically, at one-loop order one has 




The terms arising from the heavy string modes are removed by taking the low-energy limit in 
which all external momenta lie well below the energy scale set by the string tension. This limit 
picks out the regions of integration in the string diagram corresponding to particle-like graphs, 
but with different diagrammatic rules. 

Given these rules, one may formulate a purely field theoretic framework which reproduces 
them. In the case of QCD, a key ingredient is the use of a special gauge originally derived from 
the low-energy limit of tree- level string amplitudes. This is known as the Gervais-Neveu gauge 
and it is defined by the gauge-fixing lagrangian density 

£ G n = - § Tr (dpA" - i| A„ A") * . (60) 

This gauge choice simplifies the colour factors that arise in scattering amplitudes. The string 
theory origin of gauge theory amplitudes is then most closely mimicked by combining this 
gauge with the background field gauge, in which one decomposes the gauge field into a classical 
background field and a fluctuating quantum field as = A^ 1 + A^ u , and imposes the gauge- 
fixing condition D^A^^ = where is the background field covariant derivative evaluated 
in the adjoint representation of the gauge group. This hybrid gauge is well-suited for computing 
the effective action, with the quantum part describing gluons propagating around loops and the 
classical part describing gluons emerging from the loops. The leading loop momentum behaviour 
of one-particle irreducible graphs with gluons in the loops is very similar to that of graphs with 
scalar fields in the loops. 

Super symmetric decomposition. String theory also suggests an intimate relationship with 
supersymmetry. For example, at tree-level QCD is effectively supersymmetric because a multi- 
gluon tree amplitude contains no fermion loops, and so the fermions may be taken to lie in 
the adjoint representation of the gauge group. Thus pure gluon tree-amplitudes in QCD are 
identical to those in supersymmetric Yang-Mills theory. They are connected by supersymmetric 
Ward identities to amplitudes with fermions (gluinos) which drastically simplify computations. 
In supersymmetric gauge theory these identities hold to all orders of perturbation theory. 

At one-loop order and beyond, QCD is not supersymmetric. However, one can still perform a 
supersymmetric decomposition of a QCD amplitude for which the supersymmetric components 
of the amplitude obey the supersymmetric Ward identities. Consider, for example, a one-loop 
multi-gluon scattering amplitude. The contribution from a fermion propagating in the loop can 
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be decomposed into the contribution of a complex scalar field in the loop plus a contribution 
from an J\f = 1 chiral supermultiplet consisting of a complex scalar field and a Weyl fermion. 
The contribution from a gluon circulating in the loop can be decomposed into contributions of 
a complex scalar field, an N = 1 chiral supermultiplet, and an M = 4 vector supermultiplet 
comprising three complex scalar fields, four Weyl fermions and one gluon all in the adjoint 
representation of the gauge group. This decomposition assumes the use of a supersymmetry 
preserving regularization. 

The supersymmetric components have important cancellations in their leading loop momen- 
tum behaviour. For instance, the leading large loop momentum power in an n-point 1PI graph 
is reduced from \k\ n down to |/c| n ~ 2 in the Af = 1 amplitude. Such a reduction can be extended 
to any amplitude in supersymmetric gauge theory and is related to the improved ultraviolet 
behaviour of supersymmetric amplitudes. For the M = 4 amplitude, further cancellations re- 
duce the leading power behaviour all the way down to |/c| n_4 . In dimensional regularization, 
N = 4 supersymmetric loop amplitudes have a very simple analytic structure owing to their 
origins as the low-energy limits of superstring scattering amplitudes. The supersymmetric Ward 
identities in this way can be used to provide identities among the non-supersymmetric contribu- 
tions. For example, in M = 1 supersymmetric Yang-Mills theory one can deduce that fermion 
and gluon loop contributions are equal and opposite for multi-gluon amplitudes with maximal 
helicity violation. 

Scattering amplitudes in twistor space. The scattering amplitude in QCD with n incoming 
gluons of the same helicity vanishes, as does the amplitude with n — 1 incoming gluons of one 
helicity and one gluon of the opposite helicity for n > 3. The first non-vanishing amplitudes are 
the maximal helicity violating (MHV) amplitudes involving n — 2 gluons of one helicity and two 
gluons of the opposite helicity. Stripped of the momentum conservation delta-function and the 
group theory factor, the tree-level amplitude for a pair of gluons of negative helicity is given by 

n 

A(k) = e^^K\k+)I[(K\l&i)~ 1 . (61) 

i=i 

This amplitude depends only on the holomorphic (negative chirality) Weyl spinors. The full 
MHV amplitude (with the momentum conservation delta-function) is invariant under the con- 
formal group 50(4,2) = SU(2, 2) of four dimensional Minkowski space. After a Fourier trans- 
formation of the positive chirality components, the complexification SX(4, C) has an obvious 
four dimensional representation acting on the positive and negative chirality spinor products. 
This representation space is isomorphic to C 4 and is called twistor space. Its elements are called 
twistors. 

Wavefunctions and amplitudes have a known behaviour under the C x -action which rescales 
twistors, giving the projective twistor space CP 3 or MP 3 according to whether the twistors 
are complex-valued or real-valued. The Fourier transformation to twistor space yields (due to 
momentum conservation) the localization of an MHV amplitude to a genus holomorphic curve 
CP 1 of degree 1 in CP 3 (or to a real line MP 1 C MP 3 ). It is conjectured that, generally, an 
^-loop amplitude with p gluons of positive helicity and q gluons of negative helicity is supported 
on a holomorphic curve in twistor space of degree q + t — 1 and genus < I. The natural 
interpretation of this curve is as the worldsheet of a string. The perturbative gauge theory may 
then be described in terms of amplitudes arising from the couplings of gluons to a string. This 
twistor string theory is a topological string theory which gives the appropriate framework for 
understanding the twistor properties of scattering amplitudes. This framework has been used 
to analyse MHV tree diagrams and one-loop M = 4 supersymmetric amplitudes of gluons. 
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